多変量回帰・・・好きですか?
2025-06-01
古くから研究されつくされており、信頼感がある
多くの統計ソフトに入っており、行うのが簡単
解釈性が高く、分かりやすい
本当?
| Base | |
|---|---|
| RHC | 1.16 [1.01, 1.32] |
| Not CHF | 1.71 [1.29, 2.25] |
| Age | 1.03 [1.03, 1.04] |
| Num.Obs. | 5733 |
| AIC | 6882.7 |
| BIC | 7089.0 |
| Log.Lik. | -3410.366 |
| RMSE | 0.45 |
| Base | Interact | |
|---|---|---|
| RHC | 1.16 [1.01, 1.32] | 1.52 [1.02, 2.26] |
| Not CHF | 1.71 [1.29, 2.25] | 1.94 [1.40, 2.69] |
| Age | 1.03 [1.03, 1.04] | 1.03 [1.03, 1.04] |
| RHC:Not CHF | 0.74 [0.49, 1.12] | |
| Num.Obs. | 5733 | 5733 |
| AIC | 6882.7 | 6882.7 |
| BIC | 7089.0 | 7095.7 |
| Log.Lik. | -3410.366 | -3409.366 |
| RMSE | 0.45 | 0.45 |
| Base | Interact | Spline | |
|---|---|---|---|
| RHC | 1.16 [1.01, 1.32] | 1.52 [1.02, 2.26] | 1.49 [1.00, 2.23] |
| Not CHF | 1.71 [1.29, 2.25] | 1.94 [1.40, 2.69] | 1.94 [1.40, 2.70] |
| Age | 1.03 [1.03, 1.04] | 1.03 [1.03, 1.04] | |
| RHC:Not CHF | 0.74 [0.49, 1.12] | 0.75 [0.49, 1.13] | |
| rcs(age, 4)age | 1.04 [1.03, 1.05] | ||
| rcs(age, 4)age' | 0.99 [0.96, 1.01] | ||
| rcs(age, 4)age'' | 1.04 [0.89, 1.21] | ||
| Num.Obs. | 5733 | 5733 | 5733 |
| AIC | 6882.7 | 6882.7 | 6884.1 |
| BIC | 7089.0 | 7095.7 | 7110.3 |
| Log.Lik. | -3410.366 | -3409.366 | -3408.027 |
| RMSE | 0.45 | 0.45 | 0.45 |
✕
◯
Call: glm(formula = death_01 ~ swang1, family = binomial(link = "logit"),
data = rhc_prep)
Coefficients:
(Intercept) swang1RHC
0.5309 0.2248
Degrees of Freedom: 5734 Total (i.e. Null); 5733 Residual
Null Deviance: 7433
Residual Deviance: 7418 AIC: 7422
通常のアウトカム式のみで一発勝負 ここはDoubly robustは使っちゃいけない
SUTVAの原理
例えば、Matching→G computation あるいは、元々のInclusionを入れて除外したあとにIPW→アウトカム式を入れる そうすることでDoubly robust estimationとなる
基本は、ドメイン知識を入れる ただし、どのような関係性かをみるのにはAICとか、尤度比検定をしても良いかも
すべてのモデルは誤っている。しかし、そのうちのいくつかは役に立つ。
| Parameter | Odds Ratio | 95% CI | p |
|---|---|---|---|
| (Intercept) | 0.25 | 0.11, 0.56 | < .001 |
| rx [Lev] | 0.92 | 0.66, 1.28 | 0.612 |
| rx [Lev+5FU] | 0.59 | 0.42, 0.83 | 0.003 |
| age | 1.01 | 1.00, 1.02 | 0.140 |
| sex | 1.03 | 0.78, 1.36 | 0.841 |
| obstruct | 1.41 | 1.00, 2.01 | 0.053 |
| perfor | 1.00 | 0.44, 2.29 | 0.997 |
| adhere | 1.65 | 1.11, 2.47 | 0.014 |
| surg | 1.47 | 1.07, 2.00 | 0.016 |
| nodes | 1.23 | 1.17, 1.30 | < .001 |
9716.5円/Cap x 4 x 365 = 1400万円/年
医療経済的にどう考えればいい?
estimand -> 誰に?
どれくらい?